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Abstract 


For full QCD vacuum expectation values we construct an expansion in 
quark loop count and in powers of a coupling constant. The leading term in 
this expansion is the valence (quenched) approximation vacuum expectation 
value. Higher terms give corrections to the valence approximation. A test 
of the expansion is presented for moderately heavy quarks on a small lattice. 
We consider briefly an application of the expansion to quarkonium-glueball 
mixing. 
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I. INTRODUCTION 


The infinite volume, continuum limit of lattice QCD hadron masses |[T] [8|| and meson 
decay constants J4|.[3[] calculated in the valence (quenched) approximation lie not far from 
experiment. Calculations using the valence approximation, however, require significantly 
less computer time than those using full QCD. Thus, in at least some cases, the valence 
approximation can serve as a cheap, approximate substitute for full QCD. For this purpose it 
would be useful to have some way to determine quantitatively an estimate of the error arising 
from the valence approximation short of a direct comparison of the valence approximation 
with full QCD. 

A possible method for finding the valence approximation’s error is given in Ref. ||. In 
the present article, we describe an alternative form of the proposal in Ref. || which we 
believe will generally require less computer time. The expansion we describe can be applied 
to any choice of quark action but is given here only for Wilson quarks. 

In full QCD, virtual quark-antiquark pairs produced by a chromoelectric field reduce the 
field’s intensity by a factor which depends both on the field’s momentum and on its intensity. 
In the valence approximation this factor, analogous to a dielectric constant, is approximated 
by its zero-field-momentum zero-field-intensity limit 0. Our expression for the error in 
valence approximation vacuum expectation values consists of an expansion in quark-loop 
count and in powers of a coupling constant. The coupling constant expansion relies on 
ideas drawn from mean-field-improved perturbation theory [j7]. Each term in the expansion 
requires as input the quantity A/5 given by (6 rf /g 2 ) — (6 /g 2 ), where g is the gauge coupling 
constant of full QCD and rf is the dielectric constant entering the valence approximation. We 
determine A/5 analytically from mean-field-improved perturbation theory to second order 
in the coupling constant. A related calculation of A/5 without mean-field improvement is 
described in Ref. ||. The remaining work of evaluating each term in the error expansion is 
done by a Monte Carlo algorithm. 

The sum of all terms in the error expansion, in principle, gives the exact value of the 
valence approximation error for any choice of the A/5. In particular, the expansion remains 
correct independent of the accuracy of the second order perturbative expression for A/5. For 
a bad choice of A/5, however, valence approximation vacuum expectations will be far from 
their full QCD values and the error expansion will predict an error correspondly large. 

We have tested our method so far only for vacuum polarization arising from quarks with 
about 1.8 times the strange quark mass and only for a collection of Wilson loop expectation 
values. For these cases, our method of estimating the valence approximation error is signif¬ 
icantly faster than direct comparison between the valence approximation and full QCD. A 
test of the efficiency of our method for lighter quark masses and other vacuum expectation 
values we hope to return to elsewhere. 

In addition to its use as an algorithm for finding valence approximation errors, the ex¬ 
pansion we describe provides a systematic way to keep track of the quantities which need to 
be evaluated, by any method, to determine quark loop corrections to valence approximation 
vacuum expectation values. Toward the end of the present article, we present a brief, quali¬ 
tative discussion of the valence approximation and corrections to the valence approximation 
for mixing between the lightest scalar glueball and scalar quarkonium states. We will show 
that a recent attempt 0 to determine glueball-quarkonium mixing misses two of the terms 
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required for this calculation. As a consequence, we believe the calculation of Ref. 0 is not 
correct. 

In Section |IT] we introduce definitions. In Section B we construct an expansion for 
the valence approximation error. In Section we discuss the weak coupling calculation 
of the shift between the coupling constant in full QCD and in the valence approximation. 
In Section 0 we describe a trial calculations using our expansion and error estimates. In 
Section |VT] we consider the valence approximation and its corrections for glueball-quarkonium 
mixing. 


II. DEFINITIONS 

For euclidean QCD on some finite lattice, let u(x,y ) be a guage link with periodic 
boundary conditions, and let M be the coupling matrix for a single quark flavor, with 
antiperiodic boundary conditions, defined by 

M(x,y) = S xy - ~ln)u(x,y) - <Wa(1 + IfMx, y). (2.1) 

The vector jj, is a unit lattice vector in the + /j direction and the are 4x4 hermitian 
euclidean gamma-matrices. 

For nf degenerate flavors of quarks and any integrable function of the gauge fields G , 
the vacuum expectation value found after integrating out quark fields becomes 

< G > =Z~ l ( dv G det(M) nf exp(-P), 

J 6 

Z— [ dvdet(M) nf exp(—P), (2.2) 

J 6 

P= Tr[u(x 1 ,x 2 )u(x2,X3)u(x3,X4)u(x4,xi)}. 

(xi,...Xi) 

Here (3 is 6 /g 2 for bare gauge coupling constant g, v is the product of one copy of S77(3) 
Haar measure for each link variable on the lattice, and the sum in the definition of P is 
over all nearest neighbor squares (aq,.. .£ 4 ) with squares differing by a cyclic permutation 
identified. The extension of Eq. ( | 2.2| ) to vacuum expectations of products of quark and 
antiquark fields and to QCD with quarks having several different masses is not needed for 
the present discussion and will be omitted for simplicity. In the present discussion, iif can 
be either even or odd. 

The valence approximation for < G > is 

< G > v = Z- 1 j dvGexp(^-P), 

Z v = J duexp(^P). (2.3) 

Here /3 V is 6 /g v 2 with valence approximation bare gauge coupling g v . It is convenient to 
name the shift between /3 V and (3 


A (3 = (3 V - (3. 


(2.4) 
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The determination of A/3 will be discussed in Section fV[ As mentioned in Section g v and 
g may also be viewed as related by a dielectric constant rj 


g v = 9 -. (2.5) 

V 

The calculation of A (3 in Section Q is, in effect, also a calculation of g. Although g is useful 
in describing the intuitive content of the valence approximation, it will not appear directly 
in the remainder of this paper. 


III. ERROR EXPANSION 


A coupling constant and quark loop expansion can now be constructed for the difference 
between the full QCD vacuum expectation < G > and its valence approximation < G > v . 
It is convenient || to express < G > of Eq. (|2.2| ) as 

< G > = Z- 1 [ dvGexp(—P + Q) 

J 6 

Z = J du expi^P + Q), (3.1) 

Q = riftrlog(M) -^-P. 


Introducing a parameter A multiplying Q, we expand < G > in powers of A, replace A by 1, 
and get 


< G > =<G> v +J2$n(G), (3.2) 

n 

<5i(G) = < (G— < G > V ){Q- < Q >v) (3-3) 

62 (G) = < (G- < G > V )(Q— < Q > v ) 2 > v , (3.4) 

63 (G) = < (G- < G > V )(Q— < Q > v ) 3 > v - 

3 < (G— < G > V )(Q— < Q >„) >^< (Q- < Q > v ) 2 >„, (3.5) 


In a coupling constant perturbation expansion of Eq. |3.2| for the difference between < G > 
and < G >„, the quantity Q carries a single quark loop. As a consequence 6 n (G) can be 
associated with diagrams containing n internal quark loops. None of the S n (G), however, are 
simply sums of ?r-quark-loop diagrams . Each includes also, through Q, counterterms arising 
from the shift between (3 of full QCD and the screened (3 V of the valence approximation. We 
will return to this observation in Section 0 

The quantity frlog(M) in Q of Eq. [hi] we now express as a coupling constant power 
series. This series is also, formally, an expansion in powers of a gauge potential. Thus, 
as usual in gauge field theories, we use the theory’s gauge invariance to transform to a 
gauge which will tend to make the gauge potential small. For this purpose, we transform 
the gauge held u(x,y) to a euclidean lattice version of Landau gauge. The held u(x,y) 
has been transformed to lattice Landau gauge if for every lattice site x the target function 
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J2yt' r [u(x,y)] is a local maximum with respect to further gauge transformations. Generally 
there are many Gribov copies of transformations taking a particular gauge held to euclidean 
lattice Landua gauge. The vacuum expectation value of any integrable function of the gauge 
held transformed to Landau gauge is then an weighted average over the Gribov copies of 
each held. The weights depend on the particular choice of algorithm for obtaining Landau 
gauge. The gauge fixing algorithm used in our trial calculation in Section 0 is discussed in 
Ref. ]5| . Our expansion does not depend explicitly on the choice of gauge fixing algorithm 
and, therefore, on the choice of Gribov copy weighting. We have not examined to what 
degree this choice might be optimized to further minimize the gauge potential and therefore 
speed the convergence of the coupling constant expansion. 

For each hxed gauge configuration, we construct a free quark coupling matrix Mq which 
approximates the interacting coupling matrix M of Eq. |2T|. For each configuration, let z 
be the average over all lattice links of tr[u(x,y)\/ 3. Let Mq be a free coupling matrix with 
hopping constant Ko chosen to give a quark mass with agrees with the mean-field-improved |Tj] 
estimate 


— -4 = — 1 

2 Kq 2 ZK 


2 ZK, C ' 


(3,6) 


where n and k c are, respectively, the hopping constant of M and the valence approximation 
to the critical value of this hopping constant. The critical hopping constant is the smallest 
value for which the pion mass becomes zero. On the right side of Eq. |3.6| , the parameter £ 
varies with gauge configuration but k and n c do not. 

Mean-held improved perturbation theory suggests zM 0 as an approximation to M. We 
therefore express brlog(M) in the form 


trlog(M) = trlog{zM 0 [l — M 0 1 (M 0 — 2 l M))} 1 


(3.7) 


and expand to obtain the 


trlog(M) = trlog(zM 0 ) — ^ —tr{[M 0 1 (M 0 — z l M)] n }. 

n ^ 


(3.8) 


For small values of the chromoelectric potential, M 0 (M 0 — 2 l M ) is linear in the potential. 


Thus Eq. ( |3.8| ) is approximately an expansion in power of the chromoelectric potential or, 
equivalently, in powers of a coupling constant. 

To evaluate the trace in the second term of Eq. ( |3.8|) we use an ensemble of complex¬ 
valued quark fields </>. For each site x and each of the 12 combinations of spin index s 
and color index c, the corresponding (j) sc {x ) we take to be an independent complex random 
variable with absolute value 1 and probability distribution uniform on the unit circle. For 
an R element ensemble of such fields </> r , 1 < r < R, we then have 


trlog(M) = tr\og(zM 0 ) — R 1 ^ — < </> r , [M 0 1 (M 0 — z 1 M)] n (J) r > . 


(3,9) 


Here is the inner product 
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(3.10) 




on the space of complex-valued quark fields. In Eq. (|3.9| ) the quantity tr\og(zM 0 ) and the 
inverse Mj " 1 can be both be found easily in momentum space since M 0 , in momentum space, 
is block-diagonal. To multiply vectors specified in position space by Mq 1 , we transform in 
and out of momentum space by fast Fourier transforms. 

If the right side of Eq. ( |T9| ) is substituted for trlog(M) in the expression for Q in Eq. |3T| , 
Eq. |3.2| becomes a coupling constant and quark-loop expansion for corrections to valence 
approximation vacuum expectation values. The only quantity in this expansion not yet 
specified is the shift A (3. Eqs. 


and ( |).0| i are formally correct for any choice of A/3. The 
rate at which these series converge, however, will be affected by this choice. 


IV. A/3 

For valence approximation calculations of the light hadron spectrum, the up and down 
quark masses are usually taken to be equal and the corresponding k is chosen by requiring 
the pion mass to have its physical value. The valence approximation (3 V is then determined 
by setting the rho mass to its physical value. Thus, in effect, A/3 is found by requiring the 
valence approximation error in the rho mass to vanish. Since the rho mass is expected to be 
determined mainly by the low-momentum behavior of the chromoelectric held, this choice 
may be viewed as a quantitative implementation of the qualitative picture of the valence 
approximation mentioned briefly in Sect, p} A class of possible alternatives consists of choices 
of A/3 which make the error, or some approximation to the error, equal to zero for other 
quantities beside the rho mass which are determined mainly by the low-momentum behavior 
of the chromoelectric held. A convenient version of this idea for the present discussion is to hx 
A/3 by requiring zero valence approximation error for the Landau gauge gluon propagator 
at minimal nonzero momentum to hrst order in quark loops and to second order in the 
coupling constant expansion. As an additional simplification, the expectation values needed 
to determine A/3 we find using (analytic) mean-held-improved perturbation theory rather 
than by Monte Carlo. For the test case considered in Section II we obtain the same value 
of A/3 by this method as produced by the non-perturbative method of Ref. [|], 

The requirement for A/3 becomes that Si(G) of Eq. |3.3| vanishes 

< (G- < G > V )(Q- < Q > v ) > v = 0, (4.1) 

where 

G = £V(P)V(-P). (4.2) 

CfJ, 

Q = ~tr(V) - \tr(V 2 ) - AT, 

V = Mq 1 {Mq - z~ l M). 

The quantity A Cfl (p) is the Fourier transformed gauge potential 

A Cfl (p) = Y, exp (—ip ■ x)A Cfl (x). (4.3) 
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for momentum vector p with a single nonzero component p\ of 2i t/L, where L is the lattice 
period, To second order in mean-field-improved perturbation theory, the link field u(x,y), 
in M, and the plaquette action P can be approximated by 


u(x, x + p) = z{l - iA Cfl (x + )T C - ^\A clx (x + ^p)T c ] 2 } 


P = 


< trU >, 


Et 1 - + -p + -i)) 2 } 

xcfii' ^ Z Z 


(4.4) 


F ctlv (x) = A cll (x + -z>) - A Cfl (x - -u) + A cv [x - -p) - A cv {x + -p). 


where the T c are an orthnormal basis for the Lie algebra of SU(3) 

tr(T c T d ) = ^5 cd , 


and < trU > v is the valence approximation plaquette expectation value. 

The vacuum expectation values in Eqsj4~5l we evaluate by lattice weak coupling per¬ 
turbation theory to second order in the valence approximation coupling constant. This 
calculation reduces to finding the two vacuum polarization Feynman diagrams in Fig. [I]. 
These diagrams for QCD are proportional to the corresponding diagrams for a U(l) lat¬ 
tice gauge theory and thus fulfill the U( 1) theory’s Ward identities. The calculation of the 
tadpole diagram, Figure |1] (b) can thus be eliminated. We obtain 


A(3 


iV( p) 


9 rif 


4 sin 2 (7r /L) < trU >, 

1 


■[n 22 (p) — n 22 (o)] 


7 i E tr (q + p/ 2 )s(q + p)(q + p/ 2).s (g)], 

L q 


(4.5) 


where each component of q in the sum over q ranges from 7 t/L to 2tt — ti/L in steps of 2i t/L. 
The propagator S(q) and vertex T M (g) are 


S(yQ " > 1/(2 k 0 ) - i 7 fiSin(q^) - cos(q M ) ’ 

T M) = sin {%) ~ il^cos(q^). (4.6) 


The limiting value of A (3 for large L without mean-held improvement has been derived 
in Ref. f|. 


V. EXAMPLE 

As a test of our method we compared valence approximation expectations < G >„, 
Eq. (|2.3|), their one-loop errors < (G— < G > V )(Q— < Q > v ) > v , Eq. (|3.2|), and the 
corresponding full QCD expectations < G >, Eq. (0, for a lattice 10 4 with (3 V of 5.679, k 
of 0.16 and nf of 2. For a 16 3 x 32 lattice at f3 v of 5.70, Ref [|TJ gives a critical n c of 0.16940(5) 
and strange quark mass m s a in lattice units of 0.097(3). Thus k of 0.16 corresponds to a 
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quark mass about 1.8 times m s a. According to Eq. (|4.5| ), A j3 is 0.243 giving a full QCD f3 
of 5.436. For this case A/3 found by the method of Ref. H is 0.244(6). 

We used 224 independent equilibrium gauge configuration in the valence approximation 
ensemble, generated by an over-relaxed pseudo heat bath algorithm, 600 random fermion 
fields R to evaluate the trace in Eq. ( |3.9| ) and 176 weakly correlated equilibrium gauge 
conhgurations for the full QCD ensemble, generated by a red-black preconditioned hybrid 
Monte Carlo algorithm. The expansion in Eq. ( |3.9| ) was carried to order n of 10. The 
calculation of < (G— < G > V )(Q— < Q > v ) > v was not turned carefully. In particular R 
of 600 in Eq. ( |3.9| ) is much larger than its optimal vaule. The time required for the valence 
approximation and error calculation was still less than 5% of the time required by the full 
QCD calculations. 

For G we used Wilson loops IR ,... VRo consisting, respectively, of paths 1 x 1, 2 x 1, all 
rotations of steps in the directions 1, 2, 3, —1, —2, —3, all rotations of steps in the directions 
i, 2, 3, —2, —1, —3, 3 x 1, 2 x 2, 4 x 1, 5 x 1, 3 x 2, 4 x 2, and 3x3. For the 3x3 loop, lRo> 
Fig. |2] shows the predicted error < (Wo— < VRo > V ){Q~ < Q >v) >v as a function of the 
highest power n of coupling strength used in Eq. ( |3.9| ) . The error converges adequately by n 
of 7. For smaller Wilson loops, Wo ,... IR, the predicted error’s convergence as a function of 
n is comparable to or faster than the convergence shown in Fig. ||. For n of 7, Fig. [3] shows 
the relative shift of the valence approximation from full QCD (< IR > — < IR >,,)/ < IR > 
and the predicted value < (VR— < PR > V )(Q— < Q >„) >„ / < PR >. To within statistical 
uncertainties, the predicted errors agree with the true errors. 

The true errors in Fig. |3| were found from the shortest full QCD run sufficient to con¬ 
firm equilibration of < PR >,... < PRo >. Nonetheless the statistical uncertainties in 
the predicted errors are much larger than those in the true errors. If we were to run the 
error prediction algorithm long enough to obtain statistical uncertainties comparable to the 
uncertainties found by a direct comparison of full QCD and the valence approximation, it is 
possible that the computer time required by the error algorithm would become comparable 
to that for full QCD. To find the uncertainty arising from use of the valence approximation, 
however, the statistical uncertainty in the error estimate does not need to be too much 
smaller than the error estimate’s central value. Used in this way, for the set of parame¬ 
ters of the test, our algorithm takes significantly less time than the shortest possible direct 
comparison of the valence approximation and full QCD. 


VI. GLUEBALL-QUARKONIUM MIXING 


We now consider briefly the valence approximation to glueball-quarkonium mixing, cor¬ 
rections to the valence approximation to mixing which follow from Eq. ( |3.2|) and a mixing 
calculation reported in Ref. ( ||). The lowest lying glueball, according to the valence approx¬ 
imation, is stable and is expected to be a scalar. Evidence that /o(1710) is composed mainly 
of this state is given in Refs. fT0HT3ll- With quark-antiquark annihilation initially ignored, 
the lightest scalar quarkonium states are also stable. Their valence approximation masses 
and evidence for their identification with observed states are discussed in Refs. HM. 

Mixing among the ss and ( uu+dd ) /R2 scalars and the scalar glueball then occurs through 
quark-antiquark annihilation. In the valence approximation, the glueball-quarkonium mix¬ 
ing energy can be extracted [T2|JT^| from the vacuum expectation value 
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C^ff) = < g(t)s(0) >„, 


(6.1) 


operator of Ref. 

< ... > v replaced by < ... > 


where g(x,t) is the smeared zero-momentum scalar glueball operator of Ref. |T5|, with vac¬ 
uum expectation subtracted, and s(x,t) is the smeared zero-momentum scalar quarkonium 
14]]. It is convenient to define also a full QCD C{t ) by Eq. ( |6.1| ) with 
A qualitative representation of C v (t ) is provided by Figure § 
giving a typical Feynman diagram contributing to the lattice weak coupling expansion for 
C v (t). 

Assuming, for simplicity, only vacuum polarization arising from u and d quarks taken to 
have degenerate mass, the one-quark-loop correction to C v (t ) can be found from Eqs. to- 
|3.3| ). For the present discussion, we will not apply the expansion of Eq. (|3.8|) . The one- 
quark-loop error in C v (t ) becomes 


C(t)-C v (t) = T M --^T P: 

6 

T m = < [g{t)s{ 0)- < g(t)s( 0) > v ][2tr\og(M)~ < 2frlog(M) >„] >„ (6.2) 

T P = < [g(t)s(0)~ < g(t)s( 0) > V ][P~ < P >«] >v ■ 


Qualitative representations of Tm and Tp are given by typical Feynman diagrams contribut¬ 
ing to their weak coupling expansions shown in Figure 0(a) and Figure 0(b), respectively. 
Among the processes contributing to T M in Figure 0(a) are glueball-quarkonium transitions 
through common pi-pi, kaon-antikaon, and eta-eta decay channels. The quantity T P , on 
the other hand, is the counterterm, discussed in Section 0 which arises from the shift A/3 
between (3 of full QCD and the screened f3 v of the valence approximation. 

In Ref. [9J a model is proposed for mixing among the valence approximation to the 
lightest scalar glueball state and the valence approximations to the lightest scalar ss and 
(mm + dd) /\/2) states. Applied to glueball-quarkonium mixing energies, this model omits 
the leading valence approximation mixing amplitude coming from C v {t) and represented in 
Figure 0. The model includes instead only transitions through common pi-pi, kaon-antikaon 
and eta-eta decay channels. These transition do contribute to Tm- Thus the model might 
be viewed as a calculation of quark-loop corrections to the valence approximation to mixing 
if not as an evaluation of the full mixing process. The equation assumed to govern mixing 
between the valence approximation glueball and quarkonium states through intermediate 
decay channels, however, entirely ignores the counter-term Tp. No argument is offered in 
support of this omission. Ref. [0] simply assumes, without proof, a relation between full 
QCD and valence approximation propagators with no term corresponding to Tp. 

With this counter-term dropped, Eq. ( |S.2| ) gives the error in C v (t) for a version of the 
valence approximation with A/3 forced to zero. Equivalently, it is easily checked that Tp is 
the derivative of C v (t) with respect to /3„. Thus by dropping the T P term from Eq. 0 
the one-quark-loop error estimate for C v (t) is altered by approximately the increment in 
C v {t) in going from f3 v — A/3 to /3 V . For /3 V of 5.93, A/3 is known to be greater than 0.23. 
Thus a lower bound on the effect of setting A/3 to zero can be found by comparing valence 
results at (3 V of 5.70 with those at /3 V of 5.93. The data in Refs. B0 then shows that 
the one-quark-loop error estimate for glueball-quarkonium energy is changed by an amount 
equal to the entire leading valence approximation to the mixing energy obtained from C v (t). 

A cross check on the consequences for valence approximation errors of forcing A/3 to 
zero can be obtained by making this change in the error formula applied to low-lying hadron 
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masses and meson decay constants. Using the data in Refs. [|l]{|] for (3 V of 5.70 and 5.93, we 
find that with A f3 forced to zero masses and decay constants are off by as much as 45%, 
rather than by less than 10% or less than 20%, respectively, for an optimally chosen A (3. 

Thus as a calculation of errors in the valence approximation to glueball-quarkonium 
mixing energies, Ref. [I| would be expected to predict significantly larger errors than actually 
occur with an optimal choice of A (3. As we mentioned earlier, however, also missing from 
the calculation of Ref. |S[] is the leading valence approximation term which can be obtained 
from C v (t). ft appears to us that Ref. || gives neither an adequate model of the full glueball- 
quarkonium mixing process nor of the corrections to the leading valence approximation to 
this process. We believe its results are simply incorrect. 

In partial defense of Ref. |J, it might be argued that although Tp is not explicitly present 
in the relation given between valence approximation propagators and those of full QCD, Tp 
is nonetheless present implicitly. The coupling between valence approximation states and 
two-body decay channels is assumed to fall exponentially with \k\ 2 , where k is the center-of- 
mass system 3-momentum carried by one of the decay products. Perhaps this exponential 
cutoff removes from the equation of Ref. || those contributions which Tp subtracts from 
our equations. For this to hold would require a surprising coincidence since no mention is 
made in Ref. || of the need for a term like Tp in the relation between full QCD and the 
valence approximation and no attempt is made to tune the cutoff to absorb this term. The 
cutoff is introduced simply as the authors’ expectation of the behavior of coupling between 
unstable scalars and their pseudoscalar decay products. 

In addition, however, it is mentioned explicitly in Ref. |J, and supported by the tables 
giving proposed values of full QCD corrections to valence approximation masses, that the 
model of Ref. || predicts full QCD masses below valence approximation masses for those 
states which are stable in the valence approximation but unstable in full QCD. This is exactly 
the result to be expected for corrections to the valence approximation given by Eqs. (|3.2|) 
and ( [3 . 3| ) with Tp removed. For the vacuum expectation of an arbitrary G, Eqs. (PD and 


(PD s ive 


< G > — < G > v = Tm — ^-Tp, 

6 


T& = < (G- < G > v )[2tr\og(M) — < 2tr\og(M) > v \ > t 

->G 


(6.3) 


Tp = < (G- < G > V ){P- < P > v ) > r 


The contribution to the error in < G > v from the term T^ f is the incremental effect of 
the color charge screening due to a single quark loop and therefore of a decrease, by some 
amount, in the QCD effective charge. As a consequence of QCD’s asymptotic freedom 
this term shifts quantities with mass units toward smaller values in full QCD than in the 
valence approximation. From our discussion earlier it follows that the term Tp has the 
opposite effect. It shifts quantities with mass units toward larger values in full QCD than 
in the valence approximation. In fact, as might be expected from the discussion of Ref. ( TTH , 
calculations of valence approximation decay constants HR] and a recent calculation of masses 


3 show that full QCD quantities for excited states are consistently larger those of the 
valence approximation. Therefore Tp for propagators of excited states is consistently larger 
in magnitude than T^, and the model of Ref. || predicts even the wrong sign for the relation 
between masses in full QCD and in the valence approximation. It appears to us this error is 
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clear evidence that the model’s cutoff on decay momenta can not have absorbed the effect 
of the omission of Tp. 
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FIG. 1. Feynman diagrams which contribute to A/3. 
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FIG. 2. The predicted relative shift in Wilson loop W \o from its valence approximation value 
as a function of the order n in coupling constant compared to the true shift of full QCD. 
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FIG. 3. The predicted relative shift in 11 Wilson loops from their valence approximation values 
in comparison to the true shift of full QCD. 
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FIG. 4. Quarkonium-glueball mixing through quark-antiquark annihilation. 
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FIG. 5. One-quark-loop corrections to the valence approximation to quarkonium-glueball mixing. 
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